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ABSTRACT 

Altritol nucleic acids (ANAs) are a promising new 
tool in the development of artificial small interfering 
ribonucleic acids (siRNAs) for therapeutical applica- 
tions. To mimic the siRNA:messenger RNA (mRNA) 
interactions, the crystal structure of the ANA: RNA 
construct a(CCGUAAUGCC-P):r(GGCAUUACGG) 
was determined to 1.96 A resolution which 
revealed the hybrid to form an A-type helix. As this 
A-form is a major requirement in the RNAi process, 
this crystal structure confirms the potential of 
altritol-modified siRNAs. Moreover, in the ANA 
strands, a new type of intrastrand interactions 
was found between the 02' hydroxyl group of one 
residue and the sugar ring 04' atom of the next resi- 
due. These interactions were further investigated by 
quantum chemical methods. Besides hydration 
effects, these intrastrand hydrogen bonds may 
also contribute to the stability of ANA:RNA 
duplexes. 



INTRODUCTION 

In antisense technology, the antisense oligonucleotides 
have to hybridize strongly and selectively with their mes- 
senger ribonucleic acid (mRNA) complement. A variety of 
nucleic acid modifications have been synthesized for these 
purposes. Insertion of a methylene group between the ring 
oxygen atom and the anomeric carbon atom of the 
furanose ring of RNA gives altritol nucleic acid (ANA, 
Figure la). This chemical insertion has a profound effect 
on the physicochemistry and the biology of these nucleic 



acids. The nucleic acid becomes chemically and enz- 
ymatically more stable than RNA (1), while keeping 
very selective and strong hybridization properties follow- 
ing Watson-Crick rules (2). As a result of this, ANA has 
scored very well in a small interfering RNA (siRNA) 
screening assay (3). 

When carrying out the same insertion in DNA, hexitol 
nucleic acid (HNA) is obtained. Likewise, HNA-modified 
RNAs show strong siRNA effects (4). The sugar rings of 
DNA and RNA (having a furanose sugar moiety) are 
more flexible than the sugar rings of HNA and ANA 
(having a reduced pyranose sugar moiety) and hybrids 
between ANA and DNA or RNA are more stable than 
hybrids between HNA and DNA or RNA, which points 
to the importance of the presence of the OH group for 
duplex stabilization. 

Herein, we report on the structural and physiochemical 
reasons for the duplex stabilization effect of this OH 
group in an ANA:RNA hybrid, which is not observed in 
regular dsRNA's. 

MATERIALS AND METHODS 

Oligonucleotide synthesis 

Synthesis and assembly of the ANA strand were per- 
formed according to Ovaere et al. (5). At the 3'-end of 
the ANA strand, an extra phosphate group was added 
because of solid support choice. The RNA sequence was 
purchased from Eurogentec. Hybridization was accom- 
plished by titration and monitored by NMR. 

Crystallization conditions 

Crystals were formed after about 1 year by the vapour 
diffusion hanging drop method at 289 K using Crystal 
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Figure 1. (a) Representation of an ANA and an HNA nucleoside in 'C 4 chair conformation together with the RNA and DNA nucleosides for 
comparison. The numbering scheme is assigned, (b) Numbering scheme for the three ANA: RNA duplexes in the asymmetric unit. 



Screen II (Hampton Research) in a condition which 
includes 0.01 M nickel chloride hexahydrate, 1.0 M 
lithium sulphate monohydrate and 0.1 M Tris 
(tris(hydroxyl-methyl)aminomethane) pH 8.5 as buffer. 
For the crystallization drop, 1 ul of 0.5 M ANA:RNA 
decamer (duplex concentration) was combined with 1 ul 
screening condition and equilibrated against 500 ul crys- 
tallization screen. One rod-form crystal (dimensions 
0.3mm x 0.10mm x 0.06mm) was obtained and cryopro- 
tected with 40% propylene glycol. 

Data collection and processing 

Diffraction data were recorded at the Swiss Light Source 
PXIII beamline (Paul Scherrer Institute, Villigen, 
Switzerland) on a MAR225 CCD detector (100 K, wave- 
length 1.000 A, crystal to detector distance 200 mm). In 
total, 200 frames (<p increment 0.5°) covering a 100° -data 
range were collected to a resolution of 1 .96 A. Two crystals 
with similar unit cell could be found in the reflection data, 
and these data were integrated separately by MOSFLM (6). 
Nevertheless, the diffraction data originate mainly from 
one crystal and only these data could be used for structure 
solution. The data were scaled by Scala (7), and processing 
statistics is listed in Table 1. The crystal belongs to the 
orthorhombic crystal class 222 with a = 26.07, b = 42.51 
and c = 157.79 A. Systematic absences analysed by 
Pointless (7) suggested P2\2\2\ as space group. 

Structure solution and refinement 

As the crystal structure of the ANA:RNA decamer was 
already partly solved in another unit cell but suffering 
from lattice translocation defects (unpublished results), a 



decamer of this crystal structure was used as molecular 
replacement model without the extra phosphate group 
and without the last base pair at the ANA 3'-side. 
Molecular replacement by Phaser (8) resulted in a helix 
oriented almost parallel to the c-axis. As the Matthews 
coefficient (2.22 A 3 /Dalton) suggested three decamers in 
the asymmetric unit, two more duplexes were searched for. 
The three decamers are stacked on each other in a 
head-to-tail manner. 

Structure refinement was performed by Refmac (9) 
using the restraints based on the crystal structure of the 
altritol adenosine building block (Supplementary Data) 
for the ANA sugar rings and standard dictionary re- 
straints for the ANA bases and the RNA strand. The 
F 0 — F c and 2F a — F c electron density maps were carefully 
studied for any inconsistency, and bond distances, angles 
and chiral volumes were monitored. The 'Find Waters' 
option of the program Coot was used to localize 239 
water molecules (10). Disordered solvent regions are 
modelled according to Babinet's principle (11). Finally, a 
R t value of 22.44% was obtained (J? free value: 23.42%). 
Refinement statistics is listed in Table 1 . Final coordinates 
and structure factor amplitudes have been deposited with 
the Protein Data Bank (30K2) and Nucleic Acid Data 
Bank (NA0770). 

Quantum chemical calculations 

Geometry optimizations were performed at B3LYP/ 
6-311++G(2d,2p) level of theory. Initial positions of the 
C, O and P atoms were taken from the X-ray geometry. 
Selection of the starting geometries is described in the 
Results section. Starting H-atom positions were generated 
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Table 1. Data collection and refinement statistics for the crystal 
structure of the ANA:RNA decanter 
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Average mosaicity (°) 
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1302 
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239 


R { value (%) 
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53.58 
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0.017 


Bond angles (°) 


2.508 
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Values in parentheses are for the outermost shell. 



arbitrarily and later fully refined in the course of the 
geometry optimizations. The model systems carried a 
total charge of —1 and comprised two altritol moieties 
as well as a phosphate group linker. For simplicity, the 
hydroxymethyl group of the altritol unit phosphorylated 
at 3' was substituted with a methyl group. All torsion 
angles along the chain 04'-C4'-C5'-05'-P-03'-C3'-C2' 
(Supplementary Scheme SI) were fixed at the crystal 
value in the course of the geometry optimizations. The 
rest of the geometrical parameters were optimized. In 
this way, the interaction is derived in a manner consistent 
with the experimental structure (an overlay of the experi- 
mental and optimized geometries is presented in 
Figure 4c). The Gaussian09 computer code, Revision 
A.l (Gaussian, Wallingford, CT, USA) was used for 
geometry optimizations. Electronic densities derived 
from the B3LYP/6-311++G(2d,2p) wave function were 
analyzed in the frame of Bader's Atoms in Molecule 
approximation (12) with the AIM2000 code (13,14). 
Interaction energies were evaluated at several theoretical 
levels using the B3LYP/6-311++G(2d,2p) optimized 
geometries (details of the interaction energy computations 
are provided in the Supplementary Data). Relative 
energies in the presence of bulk water were obtained 
from single point calculations using the B3LYP/ 
6-31 1++G(2d,2p) optimized geometries and the polariz- 
able conductor continuum solvent model (15,16) 
(C-PCM) in the standard parameterization supported by 
the Gaussian09 program: the average surface of a tesserae 
was 0.4 A 2 and the minimum radius of the added spheres 
used to create the solvent excluded surface was 0.2 A. The 
United Atom (UA0) topological model and a scaling 
factor of 1.1 were used to define the atomic radii. 



RESULTS 

Overall helical structure 

The crystal structure of the ANA:RNA hybrid decamer 
a(CCGUAAUGCC-P):r(GGCAUUACGG) reveals a 



right-handed, anti-parallel double helix with three helices 
in the asymmetric unit (Figure 2, for numbering scheme 
see Figure lb). The electron density map reveals a discon- 
tinuous sugar-phosphate backbone throughout helical 
columns with gaps between subsequent decamers. Both 
strands are held together by standard Watson-Crick base 
pairing. In the centre of the duplex, an open space can be 
seen which classifies the ANA:RNA hybrid as an A-like 
helix. The three double helices are very similar in the 
central part and at the 3'-ends (Supplementary Figure 
SI). At the 5'-end of both the ANA and RNA strands, 
some small deviations can be observed, possibly caused by 
the absence of phosphate groups at these ends. 

Overall, the backbone torsion angles are similar within 
each strand type and typical for A-type helices 
(Supplementary Table SI). For the ANA strands, the 
standard deviations on the torsion angles are relatively 
small, whereas for the RNA strands, the standard devi- 
ations for the a- and /-angles are quite large. The y-angle 
values for the ANA strands show a remarkable pattern. 
All three CI -residues, which lack a phosphate group, have 
deviating values for the y-angle compared with the other 
residues. This can be explained by steric repulsion between 
the 05'-atom of the CI -sugar moiety and the extra 3'-end 
phosphate group of the neighbouring ANA chain. All 
bases in the ANA: RNA hybrid have an antiorientation 
around the glycosidic bond with / values being lower 
for the ANA residues than for the RNA residues. 

The helical parameters calculated by the 3DNA code 

(17) , including comparison with A- and B-type helices 

(18) , are summarized in Supplementary Tables S2-S4. 
When considering the standard deviations, the values for 
all helical parameters lie within the same ranges for the 
three different duplexes. The large X-displacement 
(—5.4 A, displacement of a base pair along its short axis, 
perpendicular to the helical axis), inclination (13.4°, the 
angle between the long axis of a base pair and a plane 
perpendicular to the helical axis), slide (—2.1°, displace- 
ment between two successive base pairs along their long 
axis) and roll (6.9°, the angle for rotation of a base pair 
around its long axis, with respect to its neighbouring base 
pair in the duplex) classify the three duplexes as belonging 
to the A-family double helices. A rather high X-disp- 
lacement and slide are observed for all duplexes 
compared with the values for A-type DNA. As exactly 
12 base pairs make one full helical turn and the twist 
angle between two successive base pairs is about 30°, the 
structure of the hybrid is more particularly organized in an 
A'-form helix (19,20). 

Supplementary Tables S3 and S4 also list base pair step 
parameters for transitions between subsequent helices. 
The twist angle for these transitions are rather low, espe- 
cially between the CD and EF duplexes where the twist is 
only 19.5°. For the X-displacement on the other hand, 
very large values are observed for the transitions of one 
decamer to another, with the largest X-displacement for 
the CD-EF transition (—9.3 A). The other helical param- 
eters for the transition steps do not show large deviations 
from the values found within the duplexes. 

For the ANA residues, the altritol sugar ring puckering 
modes are assigned to the 'C 4 chair conformation. 
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Figure 2. Side view of the asymmetric unit of the ANA:RNA decamer crystal. The right side shows the molecular structure in its contoured 2F a —F c 
electron density map (1.0 sigma level). The arrangement of the three different duplexes is assigned. Discontinuity of the helical column between 
duplexes CD and EF is clearly visible. Colour codes: nitrogen, blue; carbon, grey; oxygen, red and phosphorus, orange. Created with PyMOL (31). 



The ribose sugar rings mainly adopt the Ci'-endo con- 
formation (Supplementary Table S5). This envelope 
puckering mode is observed for all sugar rings of strand 
D. One ribose ring of strand B and three sugar rings of 
strand F reside in the C2'-exo conformation, the envelope 
form with the C2'-atom displaced from the plane to the 
opposite side of the C5'-atom. The ribose ring at the 3'-end 
of strand F adopts the CY-endo envelope conformation. 

Within each strand, the sugar puckering is closely 
related to the inter-phosphate P n — P n+l distances. The 
ANA strands have an average inter-phosphate distance 
of 5.54 A, whereas this distance is slightly higher for the 
RNA strands (5.91 A). An outlying value of 6.9 A is 
observed for the P n — P n+ \ distance between residues G9 
and G10 of RNA strand F. This large distance originates 
from the deviating CY-endo sugar puckering for the G10 
sugar ring. For the transitions between subsequent ANA 



strands, the inter-phosphate distance between the 3'-end 
phosphate group and the G2 residue is on average 6.68 A. 

Intrastrand 02 — H* • - 04' hydrogen bonds 

The size of the six-membered altritol sugar ring compared 
with the smaller flve-membered ribose can result in add- 
itional interactions between subsequent residues in ANA. 
The average distance between the 02' atom of one residue 
and the sugar ring 04' atom of the next residue is 3. 12 A 
(SD 0.24 A), where it is 3.67 A (SD 0.36 A) in the RNA 
strand. To investigate the possibility of ANA backbone 
stabilization by 02'" -- 04' interactions, the hydrogen 
atoms were added. Figure 3 visualizes positions of the 
H2' atoms using the Protonate 3D algorithm employing 
the CHARMM force field (21), consistently with the pre- 
diction made using electronic structure QM computations 



(a) 
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Figure 3. (a) Side view of a part of ANA strand E with the distances between the H2' atom of one residue to the 04' atom of the next assigned, 
(b) Side view of a part of RNA strand F. The H2' atoms point away from the 04' atom of the next residue. Colour codes according to atom type: 
carbon, grey; oxygen, red; nitrogen, blue and phosphorus, orange. Stereographic images created with PyMOL (31). 



(see later). In the ANA strands, the H2' hydrogen atoms 
mainly point towards the sugar ring 04' oxygen atoms of 
the next residue, as was found for 20 of the 30 residues. 
For the RNA strands, on the other hand, the opposite 
pattern is observed. Only six H2' hydrogen atoms out of 
30 point towards the 04'-atom of the following residue. 
These findings lead to the assumption that interactions 
between the 02' hydroxyl group of one residue with the 
04' oxygen atom of the next residue contribute to the 
stability of the ANA backbone. 

On the basis of the interatomic O • • • O distances, in the 
crystal geometry of ANA, there are three potential 
acceptor sites for the H-bond donated by the 02'- 
hydroxyl group: either the 03' or 05'-oxygens (average 
0---0 distances are 2.78 A with SD 0.14A and 3.84A 
with SD 0.42 A, respectively) of the adjacent phosphate 



group or the neighbouring 04'-oxygen (average 0 -- "0 
distance 3.12A with SD 0.24 A). 

To assess the viability of these potential H-bonding 
interactions, we carried out high-quality quantum- 
chemical computations on small model systems 
(Figure 4). For a detailed description of the computational 
models, see the Interaction energy calculations part and 
Supplementary Scheme S 1 in the Supporting Information. 
These model systems differed in the choice of the initial 
orientation of the 02'-hydroxyl group. In total, we have 
found two minima of practically identical stability (differ- 
ence of the total electronic energies is 1.4 and 0.9 kcal/mol, 
in gas-phase and in bulk water, represented by the 
COSMO continuum solvent method, respectively), with 
the 02'-hydroxyl pointing towards 04' and 03', respect- 
ively. Optimization which was started from a structure 
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Figure 4. (a and b) Two optimized geometries obtained from quantum chemical calculations at the B3LYP/6-31 1++G(2d,2p) level. H-bonding 
contact confirmed by Bader-analysis is indicated with dotted lines, (c) Overlay of the B3LYP/6-311++G(2d,2p) optimized (blue) and crystal (red) 
geometries with H-bond between the 02'-hydroxyl and 04'. (d) Model used in interaction energy evaluation. For computational details, see the 
Supplementary Data. 



with 02'-hydroxyl oriented towards 05' converged to the 
structure with 02'-H* • -04' H-bond. 

We have carried out a Bader analysis of the electronic 
density for the geometry where the 02'-hydroxyl was 
oriented towards 04' and indeed we have found a (3,-1) 
bond critical point (electron density is 0.0151, and the 
Laplacian of the electronic density is 0.0124), which un- 
equivocally indicates the existence of an 02'-H---04' 
H-bond. For comparison, we have evaluated the topo- 
logical properties of the other optimized geometry, i.e. 
the one where the 02'-hydroxyl points nearly towards 
03'. In this case, no (3,-1) bond critical point was 
found, demonstrating that no H-bond can be expected 
between 02' and 03', which is most likely because the 
03'-H-02' angle (114.7°) is too small to form a stable 
H-bond. It has been shown that H-bonds between O or 
N donors (X) and acceptors (Y) require at least an X-H-Y 
bond angle of 120° in order to contribute to the stabiliza- 
tion to a noticeable extent (22). Therefore, in light of the 
Bader-analysis and the available structural data, viability 
of an 02'-H- • -03' H-bond is much less likely than that of 
a concurrent 02'-H* • - 04' H-bond. 

To see the energetic consequences of the 02'-H- • -04' 
H-bonding interaction, we have carried out standard 
interaction energy computations at CBS(T) level (see 
Supporting Information for details) for the interaction 
of the two altritol units, by replacing the phosphate 
moiety with H atoms both at the 3'- and 5'-oxygens 
involved in the phosphodiester linkage (O-H distance for 
the added hydrogens was fixed at 0.96 A). The CBS(T) 
abbreviation stands for complete basis set electron correl- 
ation calculations corrected for higher-order electron cor- 
relation effects. Such calculations are considered of 
chemical accuracy and represent established benchmarks 
(23). The computed interaction energy, — 5.6kcal/mol, 
shows visible stabilization contribution between the two 



altritol moieties and is comparable to the gas phase 
strength of a water dimer. However, the present system 
is almost exclusively stabilized by the dispersion energy 
(roughly the difference between the MP2 and HF terms, 
that amounts to — 5.2kcal/mol using the aug-cc-pVTZ 
basis set, for details see Supplementary Table S9). The 
interaction (although being intramolecular) appears to 
be even slightly stronger than that for the ribose-zipper 
type of interactions (computed interaction energy at 
CBS(T) level is ~— 4kcal/mol), which are known to effi- 
ciently stabilize folded RNAs (24). 

The calculations presented above were specifically based 
on the experimentally refined structure of chain C of the 
ANA:RNA duplex. As an example, we have taken the 
backbone geometry from the A6, A5 dinucleotide 
segment, with an experimental 02' • • • 04' distance of 
3.15 A, i.e. which is very close to the average of the 
experimentally determined 04'-*"02' distances (3. 12 A). 
Note that despite constraining the key dihedral angles 
dictating the backbone conformation (see 'Materials and 
Methods' section), the 02' • • • 04' distance relaxed to 
3.06 A just due to optimization of valence angles and 
bond lengths. We then carried out a control computation 
of a second experimental structure (taken from strand 
C, dinucleotide step C2 and G3) with an experimental 
02'- • -04' distance of 3.57 A. This geometry can be con- 
sidered as the upper limit of the crystallographically 
observed 02'- ••04' distances. Nonetheless, after 
geometry optimization (using the same optimization 
protocol as for the other geometry)^ the optimized 
02'""04' distance decreased to 3.11 A, which is very 
close to the crystallographically found average value 
(3.12 A). Due to the similarity of the optimized geometries 
obtained from the two different starting geometries, the 
CBS(T) interaction energy results were practically identi- 
cal, see Supporting Information for further details. 
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The major groove widths, measured as the P n — P n +2 
distance minus 5.8 A to account for the van der Waals 
radii o of the phosphate groups, range from 9.5 A to 
10.8 A with an average value of 10.1 A. For the 
minor groove widths, the C4' n -C4' n+2 distances minus 
3.4 A (the van der Waals radii of two carbon atoms) 
were used as reference and have an average value of 
7.9 A with 7.3-9.0 A range. The somewhat wider 
grooves, especially the major groove, are typical for the 
A-form helix (20). 

Hydrogen bond formation between the complementary 
bases in opposite strands holds the two strands of the 
double helix together. The distances of these hydrogen 
bonds lie well within standard hydrogen bond ranges (19). 

Stacking interactions and crystal packing 

The double helical structure of the ANA:RNA decanter is 
stabilized by intra- and inter-strand base stacking inter- 
actions. A representation of the stacking interactions 
between successive base pairs in duplex AB is shown in 
Supplementary Figure S2. Inter-strand stacking is 
observed for base pair steps 1, 2, 4, 7 and 9. Duplexes 
CD and EF follow the same pattern. Base pair stacking 
interactions for the transition between different duplexes 
are mainly intra-strand (Supplementary Table S6). 

The three duplexes in the asymmetric unit stack on each 
other in a head-to-tail manner along the c-axis. By means 
of a 2-fold screw axis through the helix, six decamers stack 
on each other in the unit cell forming a helical column 



running through successive unit cells. A 2-fold screw 
axis perpendicular to the c-axis gives rise to the second 
60 base pairs column in the unit cell. The columns are 
slightly curved and are therefore not exactly parallel to 
the z-direction (Supplementary Figure S3). 

Between the ANA and RNA strands of neighbouring 
columns 16 interactions at hydrogen bond distances can 
be found (Supplementary Table S7). By these interactions, 
the side-by-side packing of the double helices is stabilized. 
Six times an RNA 02'-atom from one column is linked to 
two ANA oxygen atoms (02' and 03' or 04') of a neig- 
hbouring column. Also two hydrogen bonds linking an 
02'-atom of one column to a phosphate group oxygen 
of another column are observed. A closer view of some 
of these interactions is shown in Figure 5. 

Hydration 

For the Matthews coefficient, a value of 2.22 A 3 /Dalton 
was obtained corresponding to a solvent content of 63%. 
In total, 239 water molecules were located in the crystal 
structure resulting in an average of 7.97 water molecules 
per base pair. The program Contact (25) was used to 
analyze the interactions between nucleic acid atoms and 
ordered water molecules in a distance range of 2.3-3.5 A. 
Table 2 lists an overview of these contacts. Care has to be 
taken when considering these contacts as hydrogen bonds 
since no angle restrictions were used. 

Hydration occurs mainly at the phosphate groups, with 
some more contacts for the ANA than for the RNA 




Figure 5. Interactions between neighbouring ANA:RNA duplexes. Created with PyMOL (31). 
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Table 2. Number of contacts between ordered water molecules and 
the hydration sites of phosphate groups (OP1, OP2, 03', 05'), sugar 
rings (02', 04') and minor (02, N2, N3) and major (N4, 06, N7, 
04, N6) groove in the ANA:RNA hybrid crystal structure 



Strands 


Phosphate 


Sugar 


Minor 


Major 


Total 


All 


172 


73 


64 


71 


380 


ANA 


93 


27 


26 


31 


177 


RNA 


79 


46 


38 


40 


203 


AB 


65 


24 


26 


18 


133 


CD 


39 


21 


14 


25 


99 


111 


68 


28 


24 


28 


148 



Values for all strands, for ANA and RNA strands separately and for 
the three different duplexes in the structure are given. 



phosphate groups. The OP2 atoms pointing into the major 
groove form more contacts to water molecules compared 
with the OP1 atoms (74 versus 52), a feature that has been 
observed in DNA and RNA A-family helices (26). Water 
molecules bridging successive phosphate groups are 
present in both the ANA and RNA strands (more for 
ANA than for RNA). For the 05'-oxygen atoms, more 
interactions with water molecules are found in ANA than 
in RNA strands (15 for ANA against 8 for RNA). 

The lower hydration level for duplex CD can largely be 
ascribed to reduced hydration at the phosphate groups 
and the minor groove. More water molecules are located 
in the major groove than in the minor groove, which can 
be explained by the larger amount of possible donor/ 
acceptor sites in the major groove (27). 

The sugar moieties of the ANA strands are less 
hydrated compared with the RNA strands. A smaller 
02' • • • 04' inter-residue distance is observed for the 
ANA strands (on average 3.12 A for ANA versus 3.67 A 
for RNA) together with more 02' • • • 04' contacts at 
hydrogen bond distances (2.5-3.5 A) between successive 
ANA residues (28 for ANA and only 10 for RNA on a 
total of 30 distances). The lower hydration level of the 
ANA sugar moieties can therefore be caused by loss of 
potential hydration site due to formation of the 
above-discussed inter-residue interactions. 

In Figure 6, the solvent accessible surface area (SASA) 
of the ANA:RNA decanters is shown. The difference 
between the larger six-membered altritol sugar rings and 
the ribose sugar rings is clearly visible. To obtain the polar 
and apolar contributions to the total SASA, these were 
calculated by the GETAREA program (28) 
(Supplementary Table S8). When comparing the total 
SASA of the ANA strands (without extra phosphate 
groups) with those of the RNA strands, no large differ- 
ences are noticed. Due to the extra carbon atom in the 
altritol sugar ring, the apolar part is larger for ANA than 
for RNA strands. The polar contributions on the other 
hand are slightly smaller for ANA. 

DISCUSSION 

The altritol-modified nucleic acids were developed in 
search for oligonucleotides with an increased stability 




Figure 6. Side view of the solvent accessible surface area (SASA) of the 
three stacked ANA:RNA decamers in the asymmetric unit. The surface 
was calculated with a solvent probe radius of 1.4 A. Colour codes cor- 
respond to the underlying atom type: C, grey; O, red; P, orange and N, 
blue. Created with PyMOL (30). 



over HNAs. As a molecular-dynamics simulation 
revealed that solvation of the minor groove contributes 
to the difference in stability between a HNA:RNA and 
HNA:DNA duplex, an extra 2'-hydroxyl group in 
a-position which points towards the minor groove 
should increase the duplex stability (29). A study of the 
hybridization properties of ANA had already 
demonstrated that ANA:RNA duplexes are more stable 
than HNA:RNA duplexes (1). A crystal structure of a 
HNA:RNA decamer with a sequence analogous to that 
of the ANA:RNA decamer was solved by Maier et al. 
(30). 

The HNA:RNA decamer has no extra phosphate group 
at the 3'-end of the HNA strand, the uracil bases in the 
ANA strand are replaced by thymine bases in the HNA 
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(a) ANA 3'-end 



ANA 3'-end 



RNA 5'-end 





(b) 



HNA:RNA 




ANA:RNA 




Figure 7. (a) Superposed decamers from the asymmetric unit of the ANA:RNA crystal structure (red) and the HNA:RNA crystal structure [blue, 
Maier et al. (30)]. Side view into the minor groove (left), into the major groove (right) and a side view in-between minor and major groove (centre), 
(b) Superposed decamers from the asymmetric unit of the HNA:RNA structure (left) and the ANA:RNA structure (right). The water molecules 
bridging between successive phosphate groups are presented as spheres. RNA strands and their bridging water molecules are coloured blue. HNA 
and ANA strands with their bridging water molecules are coloured green. Created with Superpose (25) and PyMOL (31). 



strand, but in the RNA strand, the uracil bases are 
preserved. This crystal structure, with a resolution of 
2.6 A and solved in space group PA{1{1, has four 
duplexes in the asymmetric unit. In the crystal packing, 
pseudo-continuous helices containing eight decamers are 



formed. The HNA:RNA decamers fit well on the 
ANA:RNA decamers with an average positional RMS 
deviation of 0.95 A, although some deviations can be 
noticed at the 5'-ends (Figure 7a). Similar hydration 
patterns are observed in the HNA:RNA and ANA:RNA 
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structures. The OP2 atoms of the modified strands are 
much more hydrated than those of the RNA strands. 
A higher hydration level of the 05' atoms is also found 
for the modified strands compared with the RNA strands. 

By the insertion of a methylene group in RNA and 
DNA, the inter-phosphate distances are shortened affect- 
ing the hydration pattern of the modified strands. o The 
inter-phosphate distances are on average 5.54 A in 
the ANA strands and 5.55 A in the HNA strands of the 
ANA: RNA and HNA: RNA crystal structures, respect- 
ively. In the RNA strands of both crystal structures, 
these distances are on average nearly equal to 5.9 A 
(5.91 A and 5.85 A in ANA:RNA and HNA:RNA, re- 
spectively), the expected value for A-type helices. By this 
shortened inter-phosphate distances in ANA and HNA, 
many more bridging water molecules between successive 
phosphate groups are observed, almost forming a water 
spine contributing to the duplex stability (Figure 7b). The 
bridging water molecules also interact with the 05'-atoms, 
which could explain their elevated hydration level. 

For the ANA strands, the insertion of a methylene 
group introduces a new type of intrastrand interaction: 
because of the larger altritol sugar ring compared with 
the ribose ring, in the ANA strands the distance between 
the 02'-oxygen atom of one residue to the 04'-atom of the 
next is smaller compared with this distance in RNA, 
enabling 02'" - "04' interresidue interactions. Force field 
computation indicated a possibility of 02' • • • 04' 
H-bonds in the ANA strand but not in the RNA strand. 
High-level quantum chemical calculations have shown 
that an 02' • • • 04' H-bond may contribute to the stabil- 
ization of the ANA backbone. It has rather remarkable 
intrinsic interaction energy ~5.5kcal/mol with significant 
part of the stabilization coming from electron correlation 
(dispersion) contribution. Such stabilization is comparable 
with the energy of a good H-bond, whereas it is less elec- 
trostatic in nature, which may increase its relative import- 
ance in water environment. The findings reinforce the 
assumption that the higher stability of ANA:RNA over 
HNA: RNA is not only caused by extra hydration at the 
02'-atoms but also is caused by inter-residue hydrogen 
bonding interactions between the 02' atoms and 04' 
sugar ring atoms of the next residue in the ANA 
strands, whereas HNA lacks the essential 02'-atom for 
these interactions. 
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